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Abstract 


Explicit  expressions  of  the  MSA  solvation  thermodynamic  functions,  the  Gibbs  free  energy, 
energy,  enthalpy,  entropy,  apparent  molar  heat  capacities  and  p^ial  molar  volume  are  derived 
starting  &om  the  Helmholtz  free  energy.  The  thermod3mamic  consistency  of  the  MSA  solva¬ 
tion  thermodynamic  functions  are  discussed.  Some  limiting  behavious  of  the  MSA  solvation 
thermodynamic  fimctions  are  compared  with  the  Bom  theory.  The  effect  of  explicit  solvent 
structures  in  the  MSA  theory  on  the  Gibbs  and  the  Helmholtz  free  energy  is  given  special 
attention.  Model  Calculations  of  alkali  and  halide  ions  in  water  are  carried  out  and  compared 
with  experimental  data.  The  apparent  molar  heat  capacities  at  constant  volume  and  pressure, 
and  the  partial  molar  volume  are  calculated  along  an  isobaric  line  to  the  critical  region  of  a 
dipolar  liquid.  A  great  deal  of  insight  has  been  gained  on  the  behavious  of  the  solvation  ther¬ 
modynamic  properties  near  the  critical  region  of  dipolar  solvent.  It  helps  to  explain  some  of 
the  recent  experimental  observations. 


I.  Introduction 


The  theory  of  electrolyte  solution  has  been  an  active  research  subject  for  a  long  time^’^. 
Recent  efforts  are  focused  on  obtaining  structures  and  thermodynamic  functions  starting  from 
molecular  models  of  solvents  which  are  treated  explicitly®’^.  There  is  little  doubt  that  ion 
solvation  in  a  molecular  solvent  is  a  complicated  process.  Sophisticated  modehng  which  impUes 
extensive  numerical  solution  of  integral  equations'*  and  large  scale  simulations  are  needed  to  get 
accurate  structural  information  and  thermod}mamic  properties®’®’^.  However  we  believe  that 
some  simple  analyiiical  theories  can  still  capture  some  of  the  dominant  features  in  the  solvation 
process,  of  which  the  Born*  theory  is  a  beautiful  example.  It  is  the  propose  of  this  paper  to 
show  that  an  analytical  theory  known  as  the  Mean  Spherical  Approximation(MSA)  is  valuable 
to  gain  insight  on  the  solvation  thermod3rnamics. 

Indeed,  the  detailed  theoretical  and  simulation  studies®’®’^  show  that  in  spite  of  sophisticated 
solvent  re-organization  and  polarization  around  ions  the  Helmholtz  solvation  free  energy  of  an 
ion  in  a  polar  solvent  is  approximately  a  quadratic  ftmction  of  ion  charge.  This  observation 
is  strongly  supported  by  experimental  data*’*®.  One  can  understand  it  by  examing  a  well 
known  effects,  such  as,  the  dielectric  saturation**’*®’*®  near  the  ion.  For  some  systems,  the 
dielectric  saturation  of  solvent  arotmd  the  ion  may  not  be  very  important  or  may  be  cancelled 
by  an  opposing  effect,  for  example,  the  electrostriction**.  There  are  simple  relations  between 
the  internal  energy  and  the  Helmholt2  free  energy  in  some  simple  theories  such  as  that  due 
to  Bom*,  which  hold  approximately  as  shown  by  the  simulation  and  sophisticated  theories. 
For  example,  the  Helmholtz  free  energy  is  half  of  the  internal  energy.  Having  these  basic 
features  the  simple  Bom  equation  may  somehow  be  able  to  give  a  good  representation  of  the 
solvation  free  energy.  It  seems  that  there  are  indeed  ways  to  find  a  consistent  set  of  radii  to 
fit  experimental  data  of  some  solvation  thermodynamics  using  the  Bom  equation*®’*®.  On  the 
other  hand,  the  Bom  theory  is,  after  all,  an  over  simplified  theory,  in  which  the  solvent  is 
treated  as  dielectric  continuum.  The  solvent  reorientation  and  structure  breaking  that  occurs 
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when  an  ion  is  dissolved  can  not  be  properly  accounted  by  the  continuum  model.  The  MSA 
is  an  analytically  tractable  Hamiltonian  models  which  may  provide  useful  information  on  the 
solvation  of  ion  in  dipolar  solvent^^’*^  Basically  the  MSA  is  a  linear  theory  but  it  treats  the 
excluded  volume  consistent!}^.  It  has  also  the  remarkable  property  of  satisfying  the  infinite 
charge,  large  density  limits  of  Onsager^^.  The  MSA  Helmholtz  free  energy  of  charging  a  cavity 
in  a  dipolar  hard  sphere  liquid  was  calculated  by  Chan  et  al^°.  The  result  is  of  the  form  of 
the  Born  formula  with  a  modified  ion  size  dependent  on  solvent  properties-.  In  this  paper,  we 
will  obtain  other  thermodynamic  functions  based  on  the  Helmholtz  free  energy.  Indeed,  the 
MSA  gives  new  physical  interpretation  to  some  of  the  solvation  thermodynamic  fimctions,  for 
example,  the  MSA  Gibbs  free  energy  has  two  terms,  one  due  to  the  ion-dipole  interaction,  and 
another  due  to  the  solvent  re-organization,  i.e.,  the  structure  breaking  effect.  As  was  shown 
by  Garisto  and  Patey^^  that  the  MSA  Helmholtz  free  energy  is  very  close  to  the  LHNC  in  a 
wide  range  of  dipole  densities.  This  is  consistent  with  earlier  observations  that  the  MSA  gives 
reasonable  thermodynamic  properties  though  the  structure  may  not  be  very  good^ .  Recently, 
Blum  and  Fawcett  find  the  MSA  is  quite  useful  in  estimating  the  outer-sphere  contribution  to 
the  activation  parameters  for  homogeneous  electron-transfer  reactions,  and  also  in  calculating 
solvation  Gibbs  free  energy  in  various  solvents. 

Experimental  data  of  electrolyte  solutions  at  high  temperature  and  pressure  has  recently 
become  available^^’^^.  The  heat  capacities  and  partial  molar  volume  show  very  interesting 
behavior  ne<ir  the  critical  region  of  water.  There  have  been  attempts*®  to  apply  the  Bom  theory 
to  understand  experimental  observations.  It  is  found  that  the  Bom  theory  does  reasonably  well. 
In  this  paper,  the  newly  obtained  formulas  are  used  to  explain  the  experimental  observations 
using  the  MSA. 

In  the  first  section,  we  derive  expressions  for  thermodynamic  functions,  such  as,  the  Gibbs 
free  energy,  internal  energy,  entropy,  heat  capacities  from  the  Helmholtz  free  energy.  The 
thermod3rnamic  consistency  of  the  MSA  solvation  thermodynamic  functions  is  discussed.  In 
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the  second  part,  the  free  energy  and  entropy  are  calculated  and  compared  with  experimental 
data  of  alkali  and  halide  ions  in  water  at  298K.  The  dipole  moment  is  adjusted  to  fit  the  MSA 
dielectric  constant  to  the  experimental  value,  and  the  dipole  moment  so  obtained  is  considered 
to  be  the  effective  liquid  phase  dipole  moment  of  water.  We  also  calculate  the  heat  capacities 
and  partial  molar  volume  at  temperatures  and  pressures  near  the  critical  point  of  a  dipolar 
liquid. 

II.  Theory 


A.  The  solvation  Helmholtz  and  Gibbs  free  energies,  entropies  and 
chemical  potentials 


The  solvation  Gibbs  free  energy  is  the  standard  free  energy  change  of  transferring  an  ion  from 
the  gas  phase  at  1  atm  to  solvent  (mole  fraction  standard  state).  We  assume  that  the  standard 
free  energy  change  for  such  a  process  can  be  broken  down  as, 

AG®  =  AG"'®  +  AG**’®  ,  (1) 


where  AG"’®  is  the  free  energy  of  solvation  of  a  non-polar  gaseous  solute  of.  the  same  size  as 
•  * 

a  given  ion  and  AG'*’°  is  the  electrostatic  free  energy  of  solvation.  The  other  thermodynamic 
excess  functions  such  as  the  Helmholtz  free  energy,  entropy,  enthalpy  and  energy  axe  of  the 
same  form.  AG"’°  can  be  estimated  by  various  theoretical  models*®.  In  this  paper,  we  only 
consider  the  electrostatic  contribution.  The  electrostatic  part  of  the  Helmholtz  free  energy  have 
been  derived  by  constructing  a  thermodynamic  charging  process*®.  Using  the  notation  of  Blum 
and  Wei*^,  it  can  be  written  as*®’*^’*^ 


-Zi^e*(l  -  1/e®) 
<Ti  +  <^d/A® 


(2a) 
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where  Z,e  is  the  charge  of  an  ion  of  species  t,  (7i,  ad  are  the  diameters  of  ion  and  dipole 
respectively.  The  equilibrium  dielectric  constant,  of  the  MSA  is  given  by  a  single  parameter 


A® 


,0  (A®(i  -f  xyr 
16 


(26) 


A®  is  a  solution  of  following  algebraic  equation^^’^^ 

where  ^2  is  a  parameter  determined  by  the  solvent  dipole  moment 

^  =  y  wVtr 

Pd  and  p  are  the  density  and  dipole  moment  of  solvent. 


(3o) 


(36) 


The  entropy  is  given  by  the  temperature  derivative  of  the  Helmhotz  free  energy  at  constant 
volume.  To  calculate  such  a  derivative  we  have  to  calculate  (|^)v-  From  (3a),  we  get 


and 


5A®  _(e0  _  +  A®) 

(  gj,  )v  2T{2eO{l  +  Ao)2  +  Ao(A®  +  3)) 

de®  2e®(l  +  3A®) 


(4a) 


5A0  A0(1  +  AO)  ■ 

The  solvation  entropy  is  therefore  given  by 

■,aAAf\  a^Af>‘,dx\ 

A5,  =-(—^)y  = - — 

Using  the  fact  that  AAf is  a  unique  function  of  62,  from  (2)  and  (4),  we  have  explicitly 

*  .e/.o  _  -Z?e=»A®(2  +  A®)(1  -  1/£®)(2(3A®  +  l)((7,Ao  +  ad)  +  (f®  -  1)(1  +  A®)(t^) 

’  2rw^3 

where  for  later  convenience,  we  defined  a  set  of  three  W  parameters 


(46) 


(5a) 


(6) 


Wi  =  4— A®(1  +  A®)  +  2  +  3A®  , 

(^d 


=  4^J”( 

<Td 

W2  =  A®(1  +  A®)<7rf(l  -  l/e®)Viyi  , 


(7a) 

(76) 
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Ws  =  {(TiXo  +  <rd)^2e°{l  +  A«)>  +  A«(3  +  A®))  . 


(7c) 


From  the  standard  thermodynamics  relation,  the  solvation  energy  is  given  by 

Af;*'’®  =  AAf ®  +  . 


(8a) 


According  to  Garisto  and  Patey*®,  the  solvation  energy  consists  of  two  terms,  ion-dipole,  Af?*j’°, 
and  dipole-dipole,  AJ?^’®,  contributions: 


As  it  was  shown^®’*^’^^ 


AFf  •®  =  AE‘^^  -I-  AE‘Jf 


AEf/  =  2AAf  •® 


(86) 


(8c) 


Equation  (8c)  holds  precisely  for  the  continuum  solvent  model.  Accurate  theory^’^®  and  simulations®’® 
also  show  that  the  Helmholtz  free  energy  is  approximately  half  of  the  solute-solvent  contribution 
to  the  solvation  energy.  Then  from  (8a), 


and 


Explicitly  we  have 


and 


AEf’°  =  ^  Af;^’®  +  TASf’^ 


AE‘Jf  =  ^AE^’®  +  TASf’° 


,  -  l/t”) 


=  TASf  ”  -  . 


(Sd) 


(8e) 


(8/) 


(89) 


(8g)  has  been  obtained  by  Garisto  and  Patey*®  by  evaluating  the  energy  integral.  AAf’° 
is  obtained  essentially  by  calculating  the  virial  type  integral,  ie,  by  the  virial  route.  As  was 
pointed  out  by  Garisto  and  Patey  the  MSA  result  of  the  solvation  thermodynamics  is  consistent 
at  this  stage  because  one  gets  the  same  result  from  the  energy  and  the  virial  route.  As  we  see 
from  (8g)  that  AE^’°  is  always  positive,  which  corresponds  to  solvent  structure  breaking  in 
charging  an  ion  . 
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In  order  to  calculate  the  Gibbs  free  omergy,  one  starts  with  the  thermodynamic  relation 

»e/.0 


_  dpdAAj 
dpd 

Following  the  same  procedures  as  in  (5)  we  have 


^^0  A  .eio 


From  (3b),  it  can  be  shown  that 


ax\  dxo  _ 


We  then  get 


AGf •“  =  AAf^  +  TASf'°  , 


(9a) 


(96) 


(9c) 


m 


that  is, 


AG''*“  =  AEf’°  . 


(9e) 


This  is,  in  fact,  a  known  result  of  the  This  shows  again  that  the  MSA  solvation 

thermodynamic  functions  are  thermodynamically  consistent  at  this  level.  From  (8g)  and  (8f), 
we  have  exphdtly. 


O’i  +  O'd/A® 


2W3 


Comparing  with  (1)  and  (9),  we  find 


AG'*’”  ^  AAf ’®  . 


(9^) 


If  one  charges  a  cavity  in  a  dielectric  continuum  the  work  needed  is  the  same  at  constant  presure 
as  constant  volume  as  long  as  remains  unchanged,  ie, 

AG'*’"  =  AAf’®  .  (9h) 

When  a  cavity  is  charged  in  a  molecular  solvent  the  work  is  different.  The  complexity  of  the 
solvation  is  reflected  in  equation  (9f).  The  first  term  is  due  to  the  electrostatic  interaction 
between  the  ion  and  the  dipoles.  The  second  term  is  the  change  in  the  dipole-dipole  energy  due 
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to  the  ion.  As  we  pcnnted  ont  befoie,  it  is  always  positive,  which  coitesponds  to  a  structure 
breaking  effect^^.  The  Helmholtz  free  energy  is  half  of  the  first  term  in  (9f).  Only  if  the  second 
term  equals  half  of  the  first  term  the  Helmholtz  free  energy  is  not  equal  to  the  Gibbs  free 
energy. 

Consider  some  limiting  cases.  In  the  high  coupling  limit*®’^^,  where  the  dipole-dipole  inter¬ 
action  becomes  infinitely  strong,  £“  — *  oo,  from  (2a),  (9f),  we  get 

,  (10a) 

This  is  equivalent  to  the  classical  Born  result^  in  the  same  limit,  and 

A5f  °  =  0  .  (106) 

This  shows  that  strong  coupling  between  solvent  molecules  makes  the  re-organization  of  the 
solvent  around  the  ion  impossible,  and  therefore  the  entropy  change  which  is  the  measure  of 
the  microscopic  order  is  zero. 


When  the  solvent  diameter  becomes  infinitely  small,  =  0,  we  recover  the  continuum  limit 
results.  We  find  the  the  Helmholtz  free  energy  is  given  by 


^  -  l/e") 


(llo) 


which  is  precisely  the  classical  Bom  free  energy.  However  the  Gibbs  free  energy  is  found  to  be 
different  form  the  Helmholtz  free  energy  in  this  limit. 


AG* 


1,0  _ 


-Zfe^l  -  l/£°) 


O’! 


(2- 


1  -  l/£“ 


1  + 


AOO+AO) 

2e'>(l+A0)» 


■) 


and  the  entropy  is  given  by 


(116) 


.  ^,0  ^  -ZhHl  -  1A°)(2  +  A°)(3A°  +  1) 

‘  <r,r(l«(3  +  +  2£“(1  +  A»)2) 


(11c) 


This  indicates  that  the  solvent  effect,  in  this  limit,  is  still  different  from  that  of  the  continuum 
dielectric  media.  One  can  calculate  the  entropy  by  taking  derivative  of  AAf’°  in  equation  (11a) 
with  respect  to  T.  The  same  result  should  be  obtained. 
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From  Blum  and  Wei^^,  the  electrostatic  part  of  the  excess  chemical  potential  of  the  MSA 
can  be  written  as  for  an  ion-dipole  mixture  of  equal  ion  size(<T  =  <r,,  Z  =  Z,,  p  =  p,) 


elinrr\  Zido{dobo  —  d2bi<Ti/(Tj) 

(12a) 

el  +  dobr^yaf) 

(126) 

where 

^  =  4irricr^e^f{kT)  , 

'(12c) 

(12(1) 

The  b  parameters  can  be  obtained  explicitly  in  the  low  density  limit  as 

6o  =  — «/2 

(13a) 

6i  =  6|k  -f  6j/c^/2 

(136) 

62  =  6"  -1-  62«2/2 

(13c) 

where  k  is  the  Debye  screening  length 

K  —  ^ 

(Ud) 

and 

,  9d^(l  + 

'  4(l-h(7.Ao/<rd)(2-|-Ao)2  ’ 

(13c) 

•  ,  •  2  -9d2(H-AO)Vrf 

'  4.t,Ao(2  +  A»)! 

(13/) 

^  2W3<q 

(13j) 

After  some  algebra,  we  have 

2^^''  2<r.e«  ’ 

(14a) 

./  -h  Acrjf)  dib\kT 

^  Z^Pd  2irpd<T^  ■ 

(146) 

The  first  term  of  (14a)  is  simply  the  solvation  Helmholtz  free  energy. 

The  second  one  gives 

the  Debye-Huckel  limiting  law.  The  second  term  in  (14b)  is  the  chemical  potential  of  the  pure 
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dipole  liquid  and  terms  inside  the  bracket  of  the  first  term  are  related  to  the  solvation  oitropy, 
that  is, 


‘  2a.e0  ’ 

(14c) 

j^TASf-®  4bpT 

Ztpd  2irpda^j  • 

(14rf) 

The  total  solvation  Gibbs  free  energy  can  be  obtained  making  use  of  the  thermodynamic  relation 


AG*''°  =  (15) 

i 

where  AGd  is  the  Gibbs  free  energy  of  the  pure  dipole.  This  result  agrees  with  what  we  obtained 
previously,  ie,  (9b).  The  cavity  and  ideal  parts  of  the  standard  chemical  potential  are  discussed 
by  Zhou  et  al^°. 


B. 


The  apparent  molar  heat  capacity  at  constant  volume  and  pres¬ 
sure,  Cpf<,  enthalpy,  AHf'^  and  partial  molar  volume, 


To  calculate  Cff  we  have  to  take  second  temperature  derivative  of  AAf or  first  deriva¬ 
tive  of  AEf'^.  From  a  standard  thermodynamic  formula,  is  given  by 


=  ( 


dAEf''^ 

dT 


)v  = 


dAEf’°  dX° 

dxo 


(16a) 


After  considerable  algebra,  we  have 

•  Cff  =  -2A5^°  -f  -  l/e°)'((l  -I- 

+X^il  +  A«)(4(7.(1  -h  2A«)/(7d  -h  3))  -f  2Wi(l  -  1/(£°)2)(3A°  -H  1)) 

-  VF2(2W3 - +  {<Td  +  AV,)2(2A«  +  3-1-  4e'’(l  +  A°)(4A°  +  1)/A°)))  (166) 

However  Cp[f  and  are  not  very  easy  to  calculate  because  they  involve  temperature  der- 
ative  at  constant  pressure  amd  pressure  derivative  at  constant  temperature  which  can  not  be 
calculated  unless  the  equation  of  state  is  known.  Explicitly, 


^.,.0  dAHf^  dXO 

^p,i  V  QJ,  )P  ' 


(17a) 
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(176) 


,dAGt'\  _  dAGf'°  dX\ 

^  ^  dP  dX^  ^  dP 


where  is  the  enthalpy 


AHf'°  =  AGf +  rA5f'’“  =  AE''’"  +  TASf’'^ 
and  we  have  used  the  fact  that  for  the  MSA,  AGf'^  =  AEf'^.  From  (16a),  we  have 

T>e/,0  _ 


(17c) 


(17d) 


There  are  actually  a  few  ways  to  calculate  AHf'^  and  hence,  C“,-.  K  the  entropy  is  given  by 
(5a),  which  is  obtained  by  taking  the  constant  volume  temperature  derivative  of  the  Helmholtz 
free  energy,  from(8a)  and  (17c),  then 


A/f''’“  =  2AEf'^  -  AAf'°  , 
axo  =  +  ASi  )l{^)v  . 


(18a) 

(186) 


We  can  also  use  (17c)  and  calculate  the  entropy  by  taking  the  constant  pressure  derivatives  of 
Adf'^  with  respect  to  temperature. 


_  ,3^ 

It  can  also  calculated  using  the  thermodynamic  relation 


(18c) 


=  Af;*'’®  +  PV^^’° 


(18d) 


where  is  given  by  (17d).  If  the  MSA  is  a  thermodynamically  consistent  theory,  these 
three  routes  should  give  the  same  result.  Unfortunately  the  MSA  is  not  a  thermod3mamically 
consistent  theory  for  AHf'°  in  this  level,  as  we  can  see  only  some  routes  give  reasonable  results. 

As  a  first  step,  (^)p  and  (^)t  should  be  calculated.  Taking  temperature  and  pressure 
derivatives  of  both  sides  of  (3a),  keeping  pressure  and  temperature  constant  respectively,  we 


(Wo) 
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8A»  8A» 

<  8P  >’■  ’ 


(194) 


Using  the  definition  of  the  compressibilities, 


a  = 


(9  = 


(fk 

V  ’ 

(|?)t 


(20a) 

(206) 


equation  (19)  can  be  simplified 


(■gy  )P  =  (“T  +  l)(gy  )>'  . 

aA"  aA" 

(^)t  =  -l3T{-^)v  . 


(21a) 


(216) 


To  calculate  a  and  0  we  need  a  state  equation  for  the  pure  dipole  liquid  which  hais  been  given 
by  Rushbrooke  et  al^^ 

'2  16 


P.  =  +  |f^)  -  (AT  . 


(1-0^ 


18irp* 


(1  +  AO) 


The  excess  Helmholtz  free  energy  with  respect  to  the  hard  sphere  system  is  given  by 

-  -(^g(2  +  A°))^(l  +  ^T  ,  (3  + AT' 

^  54kp'  '  2  ■'■{1  +  A»)<' 


where 


n  =  pdV/NA  , 


P  =  Pdf^d  > 

^  =  e"*  ’ 

.0  6(A°  - 1) 

=  “  2  +  A»  ’ 


dm 
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d\^  (2  +  Ao)2 

Na  is  the  Avogadro  number.  (22a)  has  a  critical  point  characterized  by 


t))- 

(22a) 

(226) 

(22c) 

(22d) 

(22e) 

(22/) 

{22g) 


T:  =  .225  , 


(23a) 
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p;  =  .106 , 

(236) 

p*  =  .00704 

(23c) 

and 

P* 

‘  =.295 

kT; 

(23d) 

where  the  reduced  parameters  are  defined 

P(Ta 

p*  _ 

(24a) 

T-  I’Tci 

• 

(246) 

From  (22),  the  isothermal  and  isobaric  compressibilities  can  be  calculated,  the  details  of  which 

are  given  in  the  appendix,  and  the  results  are 

^  +  Dli-W,TQ)y  ' 

(25a) 

V 

g_  ^ 

^  +  W,T{^)y 

(256) 

where 

(N 

II 

(25c) 

poy  1  +  ^  +  ^2  _  ^3 

nBT  ’  (!-{)“ 

(25d) 

„„  2{2  +  2^-f) 

(1-0^  ’ 

(25/) 

dh%  18  aA« 

{2  +  \^Y^dT’^' 

(25^) 

The  explicit  formulas  for  the  Cp  f  and  can  be  now  written  out. 

From  (17d)  and  (21b),  we 

get 

ydfl  ^ 

(26) 

From  (18a), 

Cif  =  {2C^f  +  A5f ’«)(ar  +  1)  . 

(27) 
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Using  (17c)  and  (18c),  we  have 


AHf  "  =  AGf  "  -  . 

(28a) 

and  (17a)  gives, 

»*■’  dT"^  ^  " 

(286) 

From  (16a)  and  (21a),  we  get 

dAaf-\  dAGf’"  BX- 

(28c) 

Then  the  entropy  from  this  route  is 

i- 

and 

=  df  -  T(ar+  1)(^)t(^)p  -  T^Oif(%)p  ■ 

(28e) 

where  we  have  used  the  fact  that  C^fT  is  a  unique  function  of  b^. 

We  now 

consider  another  route,  starting  from  (18d)  and  (26),  we  have 

=  aje;*'’®  -  P^TClf  , 

(29a) 

AS*'’°  =  -P0Ctf 

(296) 

and 

off  =  C5“(aT  +  1)  -  PT0{?^)t{^)p  -  PTCi,t‘(^)p 

(29c) 

where 

( 1  srn  1  'T  1 

(30a) 

-  (p)v/r) 

/)2i,0  noy  aio 

+a({a(flj  +  D^)  +  W'<r(^)r,/-))/(^  +  Djf  -  tV.r(^)v) 

(306) 
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which  are  obtained  from  (a6)  and  (a8).  We  have  used  following  notation 


2(2  +  10^ 

(1  -  0* 


(30c) 


W. 


_  -54(gj  + 

IgyiJi'.P  (2  +  A")J  '(<riA»  +  Aj)» 


(2£“{3A“  +  1)  +  (c“  -  1)(  j^(l  +  A“)  +  A”)) 

/5ft® 

-(4£°(A“  +  1)(1  +  4A«)/A°  +  2A‘‘  +  3)(£“  -  1)(1  +  A‘‘)A“)/(2;9«W3^r)  -  i^)v/T  ,  (30d) 


All  the  solvation  thermodynamic  functions  are  now  given  in  terms  of  explicit  expressions. 
They  can  be  calculated  by  only  solving  an  algebra  equation  (3a).  The  MSA  solvation  energy, 
free  energy  and  entropy  given  in  section  A  are  thermodynamically  consistent.  However  others 
involve  second  derivatives  and  equation  of  states  are  not,  only  some  give  reasonable  results. 


III.  Model  calculations  and  discussions 


Jansoone  and  Franck^^  found  that  the  equilibrium  dielectric  constant  of  water  could  be  fitted 
fairly  well  by  the  MSA  formulo,,  (2b).  We  pick  2.8  A°  as  diameter  of  water  and  adjust  the 
dipole  moment  so  that  the  dielectric  constwt  given  by  (2b)  agrees  with  experimental  data,  ie., 
78.7  at  298K.  The  density  of  the<liquid  is  chosen  as  1.0  gram  per  cm^,  and  the  effective  dipole 
movement  so  obtained  is  about  2.22  Debye.  The  solvation  Gibbs  &ee  energy,  enthalpy  and 
entropy  of  alkalis  and  halides  is  listed  in  Table  I,  where  the  energy  is  in  kJ  and  the  entropy 
is  in  JjK.  The  second  coUumn  are  the  ion  diameters  used  in  the  current  calculations  which 
are  from  the  X-ray  electron  density  measurements^.  Comparison  between  the  e:q>erimental 
data^’^^  '^  at  298K  and  the  calculations  shows  that  agreement  for  halides  are  much  better 
tham  for  alkalis.  Corresponding  Bom  calculations  were  carried  out  with  (^^)p  =  —1-357, 
which  is  the  experimental  value  for  water  at  298K^®.  As  can  be  seen  from  (9f)  MSA  correction 
to  the  Bom  Gibbs  free  energy  is  more  than  an  effective  diameter  of  the  central  ion.  The 
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structuie  bieakmg,  ie,  solvent  le-organization,  cancels  almost  half  of  the  ion-dipole  term  giving 
a  Gibbs  free  energy  close  to  the  Helmholtz  free  energy  which  are  listed  in  the  seventh  column. 
Agreement  between  the  MSA  entropy  and  the  experimental  data  is  much  better  than  that 
of  the  Bom  approximate.  In  the  MSA  theory,  the  solvation  entropy  has  contributions  from 
the  solvent  structure  breaking  and  also  solvent  re-organization  around  the  ion.  In  the  Bom 
theory,  the  entropy  is  determined  by  the  temperature  derivative  of  the  dielectric  constant. 
It  is  probably  too  simple  to  reflect  the  complexity  of  the  solvation  process.  Table  I  shows 
that  the  MSA  solvation  entropy  may  give  better  description  of  the  stmctnre  change  associated 
with  introducing  an  ion  from  gas  phase  to  liquid.  However  the  stracture  breaking  and  solvent 
re-organization  contributions  in  the  ion-dipole  model,  may  be  over-estimated  comparing  with 
those  in  the  hydrogen  bonded  water.  So  good  agreement  with  experimental  may  result  from 
cancellation  of  errors. 

The  apparent  molar  heat  capacities  and  partial  molar  volume  are  calculated  along  an  isobaric 
line  to  the  region  near  the  critical  point  with  the  pressure  being  kept  at  139  bars.  Calculations 
were  done  for  an  ion  of  z,  =  1,  <r,  =  3.27A®  and  solvent  oi  n  =  2.2  Debye,  <Jd  =  2.5A“. 
As  we  pointed  out  before,  the  MSA  is  not  a  thermodynamically  consistent  theory  for  these 
thermodynamic  quantities,  and  only  some  routes  give  reasonable  results.  We  find  that  for  the 
system  we  have  chosen  (28e)  is  fotmd  to  give  better  results.  Since  the  equation  of  state  for  of  the 
MSA  hard  sphere  dipole  has  a  very  different  critical  point  and  isobaric  line  to  the  critical  point 
from  that  of  water  we  are  not  able  to  compaure  our  calculation  with  experimental  observation^. 
As  it  is  shown  in  Fig.  1  and  2,  some  of  the  basic  features  are  well  reproduced  by  the  MSA. 
The  and  diverge  near  the  critical  point  while  changes  slowly  with  temperature. 
Just  like  the  Bom  calculation  of  Wood  et  al^^  the  tends  to  increase  near  the  critical  point, 
but  we  do  not  see  any  evidence  that  it  is  diverging.  As  we  can  see  from  (26),  (27),  (28e)  and 
(29c)  amd  C^  f  axe  related  to  the  compressibilities  that  diverge  at  the  critical  point.  This 
hais  to  be  physical  reatson  why  and  C^f  diverge  in  the  criticad  region.  From  (16a)  is 
related  to  some  well  behaved  vairiables  at  the  critical  point.  The  MSA  theory  give  a  good  clue 
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to  the  question  Wood  et  raised  in  their  Bom  calculation  on  the  critical  behaviour  of  C^f- 
The  MSA  dipole  model  has  a  set  of  classical  critical  e3q>onents^.  It  would  be  very  interesting 
to  make  detailed  comparison  with  recent  thermod3rnamic  data^.  We  would  like  also  to  apply 
our  theory  to  other  ions,  polar  and  non-i>olar  solutes^^. 


Appendix:  Derivation  of  o;  and 


(22a)  can  be  rewritten  as 


PV  =  P^V  +  Q{h%)nRrip* 


dT 

We  have  used 

and 


dT  ^  ^  di 


From  (20a)  and  (21a),  that  is, 


V 

i'^)p  = 


we 


have 


PVa  =  PV/T  +  nfir(-Z)J«a  +  +  ^^{aT  +  l)[^)v/p-) 


It  can  be  shown  after  tedious  algebra  that 


where  W4  is  defined  by  (25c).  Solving  for  a  we  get 

pv 
nPSP 


a  = 


+W, 


^  +  W4T{§)v 


(al) 


where  P°  is  defined  by  (25d)  and  (?(5§) 

65(2  +  A«)^&o(1  +  A‘‘)^  (3  +  AO)*  16 

-  astt  +  ITW^ "  ^  ^  (TTW^^  ■  .  ^  ^ 

Taking  the  derivative  with  respect  to  temperature,  and  keeping  P  constant,  we  have 

-PV/T^r^PT(-nopl(?}L\  ■  ,  dQmM. 

PVIT  +  nRT{  Dpi^{  )p+  V^dT^^^  db^  {qj>)pIp)- 


(a3) 


(a4) 


(a5) 


(a6) 


{al) 


(a8) 
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Captions 

Table  I:  The  experimental  data  of  the  ASf  is  feom  Abraham  and  Liszi**’^,  from 

Rashin  and  Honig^®  and  AG^  calculated  based  upon  these  two.  The  energy  is  in  kj  and  the 
entropy  in  J/K. 

Figure  1:  The  MSA  and  f  as  functions  of  temperature. 

Figure  2:  The  MSA  as  a  function  of  temperature. 
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